The Linear Theory of Fluctuations Arising from Diffusional 
Mechanisms — An Attempt at a Theory of Contact Noise 

By J. M. RICHARDSON 

The spectral density is calculated for the electrical resistance when it is linearly 
coupled to a diffusing medium (particles or heat) undergoing thermally excited 
fluctuations. Specific forms of the spectral density are given for several types of 
coupling which are simple and physically reasonable. The principal objective is 
the understanding of the frequency dependence of the resistance fluctuations in 
contacts, rectifying crystals, thin films, etc. 

1. Introduction 

WHEN a direct current is passed through a granular resistance such as 
a carbon microphone or a metallic-film grid leak, or through a single 
contact, there is produced a voltage fluctuation possessing a component 
called contact noise which is differentiated from the familiar thermal noise 
component by the fact that its r.m.s. value in any frequency band is roughly 
proportional to the magnitude of the average applied voltage, and is differ- 
entiated from shot noise by the strong frequency dependence of its spectral 
density. One may regard this component of the voltage fluctuation as aris- 
ing from the spontaneous resistance fluctuations of the element in question 
if one is willing to allow the resistance to have a slight voltage dependence. 
This effect has been the subject of numerous experimental investigations, 1-8 
among which we mention in particular that of Christensen and Pearson 
on granular resistance elements. These authors (henceforth abbreviated 
as CP) arrived at an empirical formula, to be discussed presently, connect- 
ing the contact noise power per unit frequency band with the applied volt- 
age, the resistance, and the frequency for several types of granular resistance. 
Their measurements covered a range of frequency from 60 to 10,000 cps, 
and involved the variation of several other parameters, i.e., pressure. More 
recently, Wegel and Montgomery 10 have measured the noise power arising 

1 H. A. Frederick, Bell Telephone Quarterly 10, 164 (1931). 
- A. W. Hull and N. H. Williams, Pliys. Rev. 25, 173 (1925). 
3 R. Otto, Hoclifrequensteelmikund Elektroakustik 45, 187 (1935). 
•"G. VV. Barnes. Jour. Franklin Inst. 219, 100 (1935). 
6 Erwin Mever and Heinz Tliicde, E. N. T. 12, 237 (1935). 

6 F. S. Gaucher, Jour. Franklin Insl. 217, 407 (1934). Bell Sys. Tech. Jour. 13, 163 
(1934). 

7 J. Bernamont, Annates dc Pliys., 1937, 71-140. 
•M. Surdin, R. G. E., 47, 97-101 (1940). 

9 C. J. Christensen and G. L. Pearson, Bell Sys. Tech. Jour. 15, 197-223 (1936). 

10 Private communication. 
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from single contacts and have obtained results in agreement with the CP 
empirical formula down to frequencies of the order of 10 -1 — 10~ 2 cps. 

Significant theoretical work upon this problem has not been attempted 
until recently. G. G. Macfarlane 11 has advanced a theory based upon a 
non-linear mechanism containing one degree of freedom which seems to be 
in agreement with the CP law. W. Miller 12 has worked out a general theory 
of noise in crystal rectifiers. His theory is linear, contains essentially an 
infinite number of degrees of freedom, and is equivalent in many respects 
to the theory discussed in this paper; however, he has not succeeded in ob- 
taining agreement with the experimental data on crystal rectifiers (which 
satisfy approximately the CP law) for any of the specific models he used. 

The purpose of this paper is the calculation of the spectral density of the 
fluctuations of the electrical resistance when it is linearly coupled to a diffus- 
ing medium (particles or heat), or, mathematically speaking, is equal to a 
linear function of the concentration deviations of this diffusing medium. 
This diffusing medium undergoes thermally excited fluctuations and thereby 
causes fluctuations in the resistance. The motive behind this investigation 
was the understanding of the frequency dependence of contact noise dis- 
cussed in the following paragraphs, but at the present time it is apparent 
that this treatment in addition may apply to rectifying crystals, thin films, 
transistors, etc. The quantitative details of the coupling between the resist- 
ance and the diffusing medium are not considered here ; in consequence of 
which, this work can hardly pretend to give a complete explanation of con- 
tact noise. However, important results are given concerning the relation 
between the spectral density of the resistance, on one hand, and the geom- 
etry of the coupling and the dimensionality of the diffusion field on the 
other. 

Now let us consider the CP empirical formula in detail. Let R be the 
average resistance 13 of the contact (we will henceforth consider only contacts 
and will regard a granular resistance as a contact assemblage) and let Ri 
(t) be the instantaneous deviation from the average. By theorems 1-3 of 
Appendix I, we can express the m.s. value of Ri as a sum of the m.s. values 
of R\ in each frequency interval as follows: 

Rl = [ 5(«) da, (1.1) 

" G. G. Macfarlane, Proc. Pliys. Soc. 59, Pt. 3, 366-374 (1947). 

12 To be published. 

13 The resistance of a contact is composed of two parts: the "gap resistance and the 
"spreading resistance." The term "gap resistance" is self-explanatory. The "spreading 
resistance" is the resistance involved in driving the electric current through the body of 
the contact material along paths converging near the area of lowest gap resistance. The 
measured contact resistance is the sum _of these two parts. In some of the particular 
physical models considered in Section 5, R is taken to be the gap resistance necessitating 
ad hoc arguments relating gap resistance and total resistance. 
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where S{u>) is called the spectral density of R t and a> is the frequency in 
radians per second. Now in our notation the CP formula may be expressed 

5(w) = AT- 2 k b+2 /a, (1.2) 

where V is the applied d-c voltage across the contact, A is a constant de- 
pending upon the temperature and the nature of the contact, and a and b 
are constants having values of about 1.85 and 1.25 respectively. CP state 
that the constant A is equal to about 1.2 X 10 -10 in the case of a single 
carbon contact at room temperature. 

In this paper we will regard the nonvanishing of a — 2 as arising from 
a non-linear effect which should become negligible at a sufficiently low 
voltage, although this interpretation does not seem completely justified on 
the basis of the work of CP. Consequently we assume that a — > 2 as V — *■ 
in such a way that V a ~- — » 1. This is in keeping with the idea that the 
resistance fluctuations are truly spontaneous — at least for small applied 
voltages. 

Although Eq. (1.2) may represent the observations over a large range of 
frequency it must break down at very high and very low frequencies in 
order that the noise power be finite (or, in other words, in order that the 
integral (1.1) converge). 

One has several clues to be considered in looking for an underlying mecha- 
nism of the resistance fluctuations. First of all, the mechanical action of 
the thermal vibrations in the solid electrodes of the contact seems to be 
unimportant because of the following reasons: (1) there are no resonance 
peaks in S(co) at the lowest characteristic frequencies of mechanical vibra- 
tion of the contact assembly; (2) S(a>) becomes very large far below the 
lowest characteristic frequency; and (3), according to CP, Rl is strictly 
proportional to V- when the fluctuations are produced by acoustic noise 
vibrating the contact, whereas R\ is proportional to K— 2 , a ~ 1.85, when 
the fluctuations arise from the dominant mechanism existing in the macro- 
scopically unperturbed contact. One of the obvious mechanisms left is a 
diffusional mechanism. Such a mechanism does not violate any of the ob- 
servations to date and, furthermore, possesses a sufficient density of long 
relaxation times to give large contributions to 5(co) near zero frequency. 

Evidence that diffusion of atoms (or ions) can be important in modulating 
a current is provided by the "flicker effect" in which the emission of elec- 
trons from a heated cathode is caused to fluctuate by the fluctuations in 
concentration of an adsorbed layer. We might suppose that contact noise 
is a different manifestation of the basic mechanism involved in the flicker 
effect. 

In view of these considerations it seems worthwhile to investigate in a 
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general way a large class of models involving resistance fluctuations arising 
from diffusional mechanisms. In the next section we propose a general 
mathematical model embracing a class of linear diffusional mechanisms. 
In Sections 3 and 4 the consequences of the general mathematical model 
are obtained by the "Fourier" and "Smoluchowski" methods, respectively, 
these alternative methods leading to identical results. In Section 5, the 
general results are specialized to several physical cases, some of which are 
introduced only for the purpose of providing some insight into the relations 
between the possible physical mechanisms and the resultant resistance 
fluctuations, and one of which along with its refinement is a successful 14 
attempt to provide a theory of Eq. (1.2). Section 6. is a summary. 

2. The General Mathematical Model 

The physical models which we consider in this paper are concerned with 
the fluctuations of contact resistance arising from a diffusional process. We 
are consequently led to consider the following general mathematical model 
embracing a rather extensive class of the physical models as special cases: 
Let us consider the instantaneous contact resistance R{t) to be related to 
the intensity c(r, /) of some diffusing quantity as follows 15 : 

G(R(l)) = f F(r, c(r, /)) dr, (2.1) 

where r is a vector in two or three dimensional space depending on whether 
the diffusion takes place on a surface or in a volume, and dr is correspond- 
ingly a differential area or volume. The intensity c(r, /) may be either a 
concentration (in the case of diffusion of material in two or three dimensions) 
or a temperature (in the case of heat flow in three dimensions). In writing 
Eq. (2.1) we have evidently assumed that the contact resistance R(t) is 
independent of the applied voltage. Eq. (2.1) may of course allow a de- 
pendence on voltage through the quantity c; however, we will consider no 
processes involving a dependence of c on the voltage. These restrictions, 
strictly speaking, make the model applicable only in the limit of low applied 

voltages. 

Before proceeding further let us limit the treatment to the case in which 
the deviations of R and c from their average values are sufficiently small 
for higher powers of these deviations to be neglected. Let 

R(t) = R + Rd), (2.2) 

c(r, l) = c+ Cl (r, I), (2.3) 

"That is, successful in so far as agreement with the form of Eq. (1.2) is concerned. 

16 A relation more general than R{l) = SF (r, c(r, /)) dr is required as one can see from 
considering the special case of a total resistance composed of a parallel array of resistive 
elements. 
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where R and c are the average 16 values of R(t) and c(r, t) respectively. 
Evidently, R^t) = 0, and ci(r, /) = 0. Introducing the expressions (2.2) 
and (2.3) into Eq. (2.1), expanding in terms of Ci(r, I), and neglecting terms 
of the order of c\, we get 

*i(0 = / f(T) d(r, t) dr, (2.4) 



where 

c) 



M - [i% 



' dG(R) ~\ 
. dR \h--r 



The function /(r) defines the linear coupling between Ri and c\ and depends 
upon the specific physical model used. The non-linear terms neglected in 
Eq. (2.4) may be of importance under some conditions; however, we will 
not consider them here. Nevertheless, non-linear effects in the behavior of 
c\ itself are possibly important in determining the form of the power spec- 
trum of Ri(l) in the neighborhood of zero frequency. 

3. The Fourier Series Method of Solution 

In this section we consider the state of the diffusing system to be defined 
by the Fourier space-amplitudes c*(/) of Ci(r, t). The time behavior of c k {l) 
will be described by an infinite set of ordinary differential equations con- 
taining random exciting forces according to the conventional theory of 
Brownian motion. 17 This method yields the spectral density of Ri(t) directly. 

Now the diffusion process is assumed to occur in a rectangular area Ai = 
L\ X L2 or in a rectangular parallelopiped of volume A% = L\ X Li X £3 . 
In regions of the above types, if we apply periodic boundary conditions 18 , 
c\(r l) may be expanded in Fourier space-series as follows: 

Cx(r, t) = E' c k (t)e ikr (3.1) 

k 

where the components of k take the values 

ki = liziii/Li, i = 1, ••• , p, (3.2) 

in which w,- are integers and v is the number of dimensions. The prime on 
the summation indicates that the term for k = is to be omitted. This is 
required by the equivalence of the time and space averages of d (true for 
A, sufficiently large) and by the vanishing of the time average of Ci (by 
definition). 

16 The average values here may be considered as either time or ensemble averages but 
not space averages. 

17 See Wang and Uhlenbeck, Rev. Mod. Phys., 17, 323-342 (1945). 

18 If the final results are given by integrals over k- space they will be insensitive to the 
boundary conditions. 
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Before proceeding to the solution itself let us consider what it is that we 
wish to know about c k (l). Expanding the function f(r) of Eq. (2.4) in a 
Fourier space-series in the region A, , 

fir) = E V", ( 3 - 3 ) 

k 

we can write Eq. (2.4) in the form 

R l (l) = A v Z'ftc k (0 (3-4) 

k 

where ft is the conjugate of f k . 

The spectral density S(w) of Ri(t) is then 

S(«) = AlT,'C kk .(o>)ftfk', (3-5) 

kk' 

where C k k- («) is the spectral density matrix for the set c k {t) given by 
Cfcfe'(o)) = 2tt Lim - [c k (oi, r)ct-(co, t) + c k (-<a, t)c*-(-co, t)] (3.6) 

T-»0O T 

in which 

r/2 

For a full discussion of spectral densities and spectral density matrices see 
Appendix I. Consequently our objective in this section is the calculation 
of the maxtix C k k- (w) defined by Eq. (3.6). 

Now we assume that ci(r, /) satisfies the diffusion equation 

i Cl (r, t) = Dv\(r, + g(r, t) (3.8) 

dt 

where D is a constant, V 2 is the Laplacian operator in two or three dimen- 
sions, and where g(r, l) is a random source function, whose Fourier space- 
amplitudes g k (t) possess statistical properties to be discussed presently. 
The random source function g is required for exciting Ci sufficiently to main- 
tain the fluctuations given by equilibrium theory. In the case of material 
diffusion the random source function g may be discarded in favor of a ran- 
dom force term of the form -D/ X T-V-[f(c + d)], where V-f = divf, X 
is the Boltzmann constant, T is the temperature, and /is the random force; 
however, in the linear approximation these two procedures will give identical 
final results. In the case of heat flow it is understood that the diffusion 
constant is D = K/pC where K is the thermal conductivity, p the density, 
and C the specific heat. Eq. (3.8) as written is valid only for D a constant 
and C\ small. 



1 r^ 12 
c k (o>, r) = ^ / c k (l)e-"« dt. (3.7) 

lir J-t/2 
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Introducing the expansion (3.1) and the expansion 

g(r, I) £' g k (t)e ik -' (3.9) 

k 

into Eq. (3.8) we obtain the infinite set of ordinary differential equations 
-£ c k (l) = -Dk 2 c k (l) + gk (l), (3.10) 

* = | * |, 

describing the time behavior of the Fourier space-amplitudes c*(/). The 
Fourier space-amplitudes gk(l) are assumed to be random functions of / 
possessing a white (flat) spectral density matrix Ckk , independent of fre- 
quency. Multiplying Eq. (3.10) by — e~ tat , integrating with respect to time 

from — £t to +\t, and neglecting the transients at the end points of the 
r-interval, we obtain 

103 -J- UR l 

where c fc (w, t) is given by Eq. (3.7) and g k (u, t) is given by an analogous 
equation. Forming the spectral density matrices we get for the diagonal 
elements 

The matrix Gkk' can now be evaluated by the thermodynamic theory of 
fluctuations (See Appendix II). This theory gives 

c k (t)ct(t) = (3-13) 

where 

Jl if k = k' 
10 otherwise, 



\d's\ 1 \d 2 e 



\ac 2 / e =J T \dc 2 



(3.14) 



5 and e being the entropy and energy, respectively, per unit area or volume, 
T the average temperature, and x the Boltzmann constant. In the case 
where c is the concentration of particles whose configurational energy is 
constant, s" = x/c. If c be the temperature T then s" = C/T- where C 
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is the heat capacity per unit area or volume. Now by a general theorem 
concerning spectral density matrices (see Appendix I) we have 



CkiOctiO = / C k k'(co) dot, 
Jo 



giving finally by combination with (3.12) and (3.13), 



, 2 



Gkk' = - -. — ti ^kk' , (3.15) 

7T A v S 



and 



1.2 



Cfcft'(oj) = --r—r, „ , n2 , 4 - (3.16) 

7T A,S CC* -]- V k* 



The spectral density S(u>) of i?i(/) then becomes 



,2 \ f |2 



^) = ^Z^L 4 . (3.17) 

If we are concerned with frequencies greater than a characteristic fre- 
quency 

wo = 47r 2 ZyZ, 2 (3.18) 

where Z, is the smallest of Li , i = 1, • • • , y, then the summation in (3.17) 
may be replaced by an integration giving 

c , * r+ i ,-ixD f \f(k)\ 2 k 2 dk 

S(») = 2 tt -p, J w2 + jD2 ^ 4 (3.19) 

where 

^Z7r; Jx r 

The integration in Eq. (3.19) is carried out over the entire ^-dimensional 
fe-space. If the range of the function /(r) is sufficiently small compared with 
the region A v , or if we let A t become indefinitely large, then the integration 
in Eq. (3.20) may be extended to all of ^-dimensional r-space. 

It is perhaps revealing to rephrase Eqs. (3.17) and (3.19) in terms of 
distributions of relaxation times. In the theory of dielectrics we speak of 
the real part of the dielectric constant being equal to a series of terms 
summed over a distribution of relaxation times: 2< a»Ti/(l + Tj&> 2 ), if the 

distribution is discrete, or / c(t)t(Zt/(1 + tV), if the distribution is con- 
Jo 

tinuous. In the above, a, is the weight for the relaxation time n , and, in 
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the case of a continuous distribution, a(r)dr is the weight for the relaxation 
times in the range dr containing r. In these terms Eq. (3.17) becomes 

5(«)=i;-^f- 2 , (3.17a) 

k 1 + TfeCO 



where 



Eq. (3.19) becomes 



a* = - i |A | 2 . (3.17b) 

IT S 



where 



SM = f £t± (3.19a) 

JO 1+ T-OP 

«W - p^T^i / l/C/V^) I 2 «- (3.19b) 



in which Z is the unit vector in the direction of k, dft, is the differential 
"solid" angle in the p-dimensional /j-space, and the integration is over the 

total solid angle (2ir in 2 dimensions, or 4ir, in 3). 

It is of interest to calculate the self-covariance Ri(t)Ri(t + u). In Appendix 
I, it is shown that the self-covariance above is related to the spectral density 
S(oi) as follows: 

Ri(i)Ri(l + u) = [ S((a) cos tua da. (3.21) 

Jo 

Using S(u) in the form (3.17), Eq. (3.21) gives 

RiiORrO + u) - xA,/*" ■ £ I/* !*«■*■*, ( 3 - 22 ) 

ft 

« > 0; 

whereas with 5(<o) in the form (3.19) we get 

MMt + «) = (2t)' xA" / I fik) | 2 <T° ufc2 <fc. (3.23) 

The method of the next section yields the self-covariance directly. 

4. Smoluchowski Method of Solution 

We call the procedure employed in this section the "Smoluchowski 
method" because it is based on an equation very closely analogous to the 
well-known Smoluchowski equation forming the basis of the theory of 



126 BELL SYSTEM TECHNICAL JOURNAL 

Markoff processes. 19 We set out directly to calculate the self-covariance for 
Ri{t) which, by Eq. (2.4), is given by 

Rx(t)Ri(t + «) = fjf(r , )f(r)c 1 (r't)c 1 (r , ,t + u) dr' dr. (4.1) 

Thus the problem is now reduced to the calculation of Ci(r', t)ci{r, t + «)• 
The quantity Ci(r', /)ci(r, I + u) is calculated in two steps. First we find 

Ci(r, / + u) , the average value of c x at the point r at the time t 4- u 
with the restriction that d is known at every point r' with certainty to be 
d(r', /) at the time / (assuming, of course, that u > 0). Then we find that 
the required self-covariance for d is given by multiplying the above 

Ci(r, t + u) by C\{r' , t) and averaging over-all values of Ci(r, t) at 

time /; thus: 



ciif, (Mr, t + uf +u) = ci(r', /)ci(r, t + «). (4.2) 



.(t+u) 



Now we assume that c x (jr, I + u) is related to c(r', t) by an integral 

equation, analogous to the Smoluchowski equation, as follows: 

('+«) r 

ci(r, t + u) = J P (\ r - r' |, u)d(r', t) dr'. (4.3) 

In the case that c represents a concentration as in the diffusion of particles, 
p(| r — r' |, u) dr is the probability that a particle be in the v-dimensional 
volume element dr at time / + u when it is known with certainty to be at 
r' a time /. Now the number of particles in dr' at f at time / is evidently 
[c + Ci(r', /)] dr'; consequently, the probable number of particles in dr at 
time t + u which were in dr' at time / is p(| r — r' \, u)[c 4- Ci(r', /)] dr dr'. 
Integration over dr' gives the total probable number 

(c + Ci(r, / + m) ) <fr of particles in dr equal to 

\\p{\ r ~ r ' I) M )[* + c i( r '» 01 <fr' ) ^ r which reduces to 

( 2 + / p{\ r — f |, u)ci(r', /) t/r' j dr. Division by </r and subtraction of c 

from both sides of the equality leads directly to Eq. (4.3). For the case of 
heat flow in crystal lattices the above picture can be used approximately if 
one uses the concept of phonons. 20 For a diffusional process p(| r — r' |, u) 
is the normalized singular solution of the diffusion equation 21 ; thus 

p {\r -r'\,u) = ^ exp [- 1 r - r' \*/4Du] (4.4) 

"Loccit 

2 "J. Weigle, Experientia, 1, 99-103 (1945). 

21 Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943). 
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where v, as previously defined, is the number of dimensions of the region in 
which the process occurs. 

Combining Eqs. (4.2) and (4.3) we get 

ft(r', iMr, l + «) = f W, Ociir", t) p(\ r - r" |, u) dr". (4.5) 

Now using the fact that 

(O . 



ci(r', tMr", i) = ctf, t) Cl (r", i) (4.6) 

and using the relation 

*(r', /)*(r", t) = i a(r' - r") (4.7) 

proved in Appendix II, Eq. (4.5) reduces to 



ci(r', t) Cl {r, l + u) = p P (\ r - r'\, «). (4.8) 

Introducing the expression (4.8) into Eq. (4.1) we obtain at once the desired 
result 

RMRiit + u) = \ fff(r)f(r%(\ r - r' \, u) dr dr' 

(4.9) 

" A4^P If f(r)f(r ' } eX P [ -' r - r ' \ /ADu] * dT '- 

For the sake of comparison with Eq. (3.23) it is necessary to write (4.9) 
in terms of the Fourier space-transforms of the pertinent quantities. We 
write 



where 



/(r) = ff(k)e k - r dk 

- m = ih-y\ f{r)eik ' rdr - 



(2*0 

Also, we write 

p(| r-r'\, u) = (4 ^ m) „, exp [-\r - r' f/Wu] 

= tJ-t f exp l-Dnk- + ik-(r - r')] dk. 
(2ir)' J 

After introduction of these expressions into (4.9) a short calculation yields 
the result 

RdORiU + «) = (2tt) v * J \f(k) \*e- Duk ' dk (4.10) 
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which is identical with Eq. (3.23) (provided that we let A v — > oo in the lat- 
ter). Thus the methods of approach used in Section 3 and in this Section 
are completely equivalent. 

5. Special Physical Models 

In the previous two Sections we have developed by two different methods 
the consequences of the general mathematical model discussed in Section 2. 
Here we apply the general results to some special physical cases. In this 
task we will be principally concerned with finding the form of the function 
fir) and establishing the number of dimensions v of the diffusion field. The 
main objective here is to provide some orientation on what mechanisms are 
or are not reasonable and to find at least one mechanism leading to the 
observed spectral density (inversely proportional to the frequency). 

a. A General Class of Models. Here we consider all at once mechanisms 
which can be adequately represented by having /(r) a v-dimensional Gaussian 
function of the form 

where Aj /2 is the "width" of the function measured along the i-th coordinate 
Xi. This form of /(r) can represent approximately several types of localiza- 
tion of the coupling between i?i and c x , as will be seen in the special examples 
later. Now if we work with A „ = oo , we will then have to consider the Fourier 
space- transform of f(r), which is readily shown to be 



(5.2) 



fW=tk'I f{r)e ~ ik '' dr 
= ft V2ir-<f A < fcf/! 

Inserting this result into Eq. (3.19) we obtain immediately 

* {a) - ^ \h=\ Tit)]- J + D> k* (5.3) 

Inserting this expression (5.2) into Eq. (3.23) gives 

ft 'W + "''? Lja)i^ " >0 (14) 
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In order that (5.3) give the observed S(co) a 1/co as a result, the integral 

nust reduce to something proportional to / k dk/ico 2 + D 2 k A ). It is clearly 

impossible that any choice of v and any set of Aj can achieve this result. 
Furthermore the self-covariance R x {i)R x {t + n) corresponding to the ob- 
served S{u>) should depend explicitly on the way S(a>) deviates from l/o> 
as w goes to zero. The expression (5.4) is finite for all u > and does not 
depend upon any cut-off phenomena in S(co) at low frequencies. Therefore 
we can exclude any physical mechanisms belonging to the class considered 
here. However, since several mechanisms that have been proposed do fall 
into this class, we consider them below: 

(t) ScliiJJ's Mechanism. SchifF considered tentatively that the fluctua- 
tions in contact resistance may be due to the variation in concentration of 
diffusing ions (atoms, or molecules) in a high resistance region bounded by 
parallel planes of very small separation. Schiff arrived at a noise spectrum 
proportional to 1/w but at the expense of disagreeing in a fundamental way 
with thermodynamic fluctuation theory. Here we will show what the cor- 
rect consequences of this mechanism are. 

Consider that the high resistance region is bounded on either side by planes 
parallel to the (x\ , :v 2 )-plane and that the thickness in the x 3 -direction is 
very small. Now this is obviously a case of the general model just considered 
in which we take v = 3 and 

Ai l2 »£>«,] 

(5.5) 
A 3 «£>«,) 

where \/u is of the order of magnitude of the frequencies of interest. It is 
then a matter of algebraic manipulation to show that 



, v x£ b\ b\ 1' 3 r kjdh 

^^ s" ' 27r 3 Ai /2 Ai /2 ' Jo » a + J *{ 



X bl b\ b\ 1 



(5.6) 



D^s" 2 5/ VAi /2 A 2 /2 " ,/2 



and 



W)W + -) = ^ ■ (ijjpg/s • Jb - > o. (") 

Thus we see that Shiff's mechanism leads to a noise spectrum proportional 
to l/co 1/2 , not 1/ft). The explanation of the singularity of the self-covariance 
(5.7) at u = lies in the inequalities (5.5). 

22 L. I. Schiff, BuShips Contract NObs-34144, "Tech. Rpt. #3", (1946). Before the 
publishing of this paper, Schiff informed the writer that he has discarded this mechanism. 
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The above treatment could just as well be applied to the case in which 
the diffusing quantity is heat instead of ions. 

(it) Resistance of a Localized Contact Disturbed by a Diffusing Surface 
Layer. Here we consider the case of two conductors covered with diffusing 
surface layers. It is supposed that the conduction from one conductor to 
the other is distributed Gaussianly with a width A 1/2 . Finally, it is supposed 
that the conductivity through the above area varies with the surface con- 
centration of the surface layer in that region. This situation is well repre- 
sented by the above general model by taking v = 2, Ai = A 2 = A, and 
b Y = b 2 = b. 

The self-covariance is readily calculated with the result 

b* 



*(04« + «)=£■ MA + Du) (5.8) 



The corresponding spectral density is 

, , 1 x& 4 r COS Ucadu 1 X& 4 T ( KirC\r-( a /n\ 



+ sin (o>A/D) ( | - Si(a)A/D) 



)] 



(5.8a) 



where Ci(x) and Si(x) are the cosine and sine integrals 23 respectively. When 
w « D/A 

S ( u ) c* - -L X* log ^A/D), (5.8b) 

Z7T" 5 L) 



and when co ^> D/A 



8 = 0.5772, 



*«~2^- <W 



Thus we see that this case does not lead to the experimental form of the 
spectral density. It must be remarked that here S(a)) is very sensitive to 
the form of the self-covariance for small u. 

b. Contact between Relatively Large Areas of Rough Surfaces Covered with 
Diffusing Surface Layers. We consider this case in detail since it leads to 
results in agreement with experiment. Furthermore, the more detailed 
consideration of this case will illustrate more fully the use of the general 
mathematical model, which may be of use in studying other diffusional 
mechanisms should they be postulated at some future time. This mechanism 
does not fall into the class just considered. 

» See Jahnke and Emde, "Tables of Functions," p. 3, Dover (1943). 
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Suppose that the contact in an idealized form consists of two rough sur- 
faces close together. Let positions on the surfaces be measured with respect 
to a plane between the surfaces, which we will call the mid-plane. Let the 
coordinate system be oriented so that the X\ and .r 2 axes lie in the mid- 
plane. Furthermore let the region in the mid-plane corresponding to close 
proximity of the rough surfaces be a rectangular area A 2 = IiX L 2 . Now, 
for convenience, we describe positions on the mid-plane by a two dimensional 
vector r = (xj , .r 2 ), and henceforth it will be understood that all vector 
expressions refer to this two-dimensional space. Let the distance between 
the upper and lower surfaces at r be denoted by h(r). The geometry of the 
above model is illustrated in Fig. 1. 

Now suppose that each surface is covered by a diffusing absorbed layer, 
such that the sum of the concentrations on both surfaces is c(r, I) at the 
time / in the neighborhood of r. Now consider the conduction of current 
between the surfaces. Let us assume that the conductance per unit area 
(of mid-plane) is a function of the separation A of the surfaces and the total 
concentration c of absorbate near the point in question, i.e., F(h, c). The 
total conductance will be the sum of the conductances through each element 
of area: hence, the instantaneous resistance R(l) at time I will be given by 

l/R(l) = f F(h(r), c{r, /)) dr (5.9) 

where dr is the differential area on the mid-plane and the integration extends 
over the rectangle A 2 = U X L 2 . Behind the above statements lies the 
tacit assumption that the radii of curvature of the rough surfaces are gen- 
erally considerably larger than the values of h. However, we will not explic- 
itly concern ourselves with this implied restriction. 

At this point it is expedient to imagine that we have an ensemble of con- 
tacts identical in all respects except for different variations of the separa- 
tion h{r). If we have any function of //, \f/(h) say, which we wish to average 
with respect to the variations of //, we simply average the function over the 
above ensemble giving a result which we denote by }f/(li) ' . 

Now let us write 



and, as before, 



h(r) = h u) + } h {r), (5.10) 

R(i) = R + *,(/), ] 

(5.11) 
c(r,l) = g + <r,(r,/). 



We assume that the ensemble average h (e) and the time averages R and c 
are constants independent of r and /. Let us also assume that the integrals 
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of fa(r) and c x (r, l) over A 2 vanish. Inserting (5.10) and (5.11) into (5 9) and 
expanding, we get 

l/R - *,(/)//? + ••■ = A t F{h u \e) 



+ 



(©■jL^.^+KS)'^* 



+ 



(5.12) 



where the super zero on the derivatives indicates that they are evaluated at 
h = J5 (e> and c = c. In accordance with previous approximations in this 
memorandum we neglect 24 terms of the order of c\ and Rl . We also neglect 
terms of the order of hi . After taking the time average of (5.12) and sub- 
tracting the result from (5.12) we get 



Ri(0 = / /(r)ci(r, dr, 

/(r) = aR 2 h(r), 



01 = ~wJ 



(5.13) 



Thus we now have a special case of our general mathematical model for the 
number of dimensions v = 2, provided that we assume that the total con- 
centration c on both of the rough surfaces fluctuates in the same manner as 
the concentration of a single adsorbed layer confined to a plane rectangular 
surface. The spectral density S(a) of &(/) is then given by Eq. (3.17) which 
we repeat here 



S(a>) = - 



2 X AoD 



2' * 2 IM 2 

k tf + DW 



(5.14) 



where k is a two-dimensional vector whose components take the values 
ki = 2*nilU , Hi = 0, ± 1, ± 2, • • • , and where f k are the Fourier space- 
amplitudes of /(r) given by 



h^AT 1 [ f(r)e- ikr dr. 

Ja 2 



It may be appropriate at this point to consider the quantity s" in detail 
for this particular case. If the energy e per unit area is independent of c, 

-2 

we have s" = — — B evaluated at c = c where s is here the entropy of the 
dc- 

absorbate per unit area. For the sake of illustration let us consider a single 
layer of absorbate in which the molecules are non-interacting. If c, the sur- 

24 We neglect these terms not because they are small compared with c\ or h\C\ but, 
because they are non-fluctuating (in time), are hence to be compared with l/R. 
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face concentration of the absorbate, be measured in molecules (atoms, or 
ions) per unit area, then, for the ideal system above, it follows that s = — x 
c log c and finally that s" = x/c. However, in the mechanism discussed in 
this part we have a compound system consisting of two separate layers on 
the upper and lower surfaces respectively. Nevertheless, a detailed analysis 
reveals that with c equal to the sum of the concentrations of both layers 
we still have s = xA even though s itself is no longer given by an expres- 
sion the same as that above. In conclusion let us consider the factor x/s" 
in Eq. (5.14). This factor is under the above idealization simply equal to 
c. That is, the spectral density S(w) is directly proportional to the average 
concentration of absorbed molecules, meaning simply that each molecule 
makes its contribution to the resistance fluctuations independently of the 
others. Of course, in any real system this will not be quite true; however, 
the existence of interactions will be manifested only by making x/s" not 
equal to c in Eq. (5.14). 

The results quoted thus far apply to a system with a given h(r). Now we 
shall average the right-hand side of Eq. (5.14) over the ensemble of varia- 
tions of h(r), it being supposed that S(co) itself on the left-hand side will 
be negligibly affected by this operation. This amounts to replacing \f k | 2 

(e) 

by I/* I 2 • We then have 

where //* are the Fourier space-amplitudes of //i(r). 

We now consider more closely the problem of calculating | hk | 2 ". We want 
to assume that h\{y) is a more or less random function of r. If //i(r) were a 
random function of r in the same way that the thermal noise voltage is a 
random function of the time /, then | li k | 2 would be a constant independent 
of k and the self-covariance hi(r)hi(r') would vanish for r ?± f. This 
clearly cannot be so, since the function lh{r) with such statistical properties 
would represent a highly discontinuous type of surface incapable of physical 
existence. We then fall back upon the more reasonable assumption that the 
gradient of hi possesses statistical properties of the above type. This notion 
is precisely formulated by means of the following equations: 

V//i(r) = p(r) (5.16) 

where 

f p(r) dr = 0, (5.17) 

and 

p(fW) W - fti«(r - r'). (5.18) 



fkf M = aW \h k \* (e) (5.15) 



•w 



134 BELL SYSTEM TECHNICAL JOURNAL 

In Eq. (5.18) /3 is a parameter (with the dimensions of area) characterizing 
the amplitude of the surface roughness, and 1 is the unit tensor in two di- 
mensions. Expressing (5.16) in terms of the Fourier space-amplitudes h k 
and pk of hi and p respectively, we have 

-ikh k = p k , (5.19) 

giving finally 

PEF W = **3S W • V* 4 ( 5 - 2 °) 

Expressing (5.18) in terms of Fourier space-amplitudes we get 

ftff, W = l3A7 ] n kk ., (5.21) 

which, when inserted into (5.20) gives the following desired result: 

rKp w = mi 1 *~ 2 - ( 5 - 22 ) 

Now replacing |/ fc | 2 by] f k | 2 ' in Eq. (5.14) and substituting the ex- 
pression (5.22) with the use of Eq. (5.15), we obtain 

7T S ft CO" + L^* 4 

If the frequencies of interest are larger than a certain characteristic frequency 
coo = 4tt D/U where L is the smaller of Li and L 2 , the summation in 
(5.23) may be replaced by an integration giving finally 



S(«) = xD/rs w -a t 0&Ar [ - s 

Jo CO 



£d£ 



„2 _f_ ^/fe* 

= x/47r/'-/3o: 2 i? 4 >l2-l/co 



(5.24) 



This result is in agreement with experiment in most respects. The dependence 
on frequency is, of course, that experimentally observed by all investigators. 
The non-dependence on the voltage applied across the contact is implied 
by the basic assumptions common to all of the mechanisms considered here, 
and is in approximate agreement with the results of Christensen and Pearson 
(see Eq. (1.2)). For our result to agree with the results of CP as regards the 
dependence on the average resistance 25 /?, the factor c 2 K* A 2 must be pro- 
portional to R 2+b where b ~ 1.25. These authors also imply that some of 

"It must be remembered that the resistance R in the CP formula is the total contact 
resistance equal to sum of the gap resistance and the spreading resistance, whereas the R 
in our theory evidently should be considered the gap resistance. For the purposes of com- 
parison we make the ad hoc assumption that the gap resistance is proportional to the total 
contact resistance. 
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the parameters necessary to complete the description of a contact between 
given substances at a given temperature show up implicitly only through 
R. According to our theory the factor or R A At does not depend in any unique 
way upon R; it matters by what means R is varied. If the resistance R is 
changed by altering the contact area A 2 , keeping other parameters fixed, 
we would find that RAo is constant so that the factor in question would be 
proportional to R 3 , that is, 6=1. However, if R is changed by varying the 
contact pressure, the effect would show up through the factor a 2 , (0 also, 
to some extent, perhaps) and, since one would expect a to increase somewhat 
with pressure whereas R decreases with pressure, the factor of interest would 
probably depend upon some power of R between 3 and 4, that is, 1 < b < 2. 
The theory formulated here suffers from the difficulty that the integral 
of the power spectrum with respect to frequency is logarithmically divergent 
at and °°, that is 

I S(oi) do> = / do: oj = log (012/(^1) — * °° as oji — > and «2— > °° . 

The divergence at » does not bother us as much as the divergence at 
since, with only a divergence at «>, the self-covariance i? 1 (/)i? 1 (^ + u) exists 
for all non- vanishing values of u; whereas, with a singularity at 0, the self- 
covariance does not exist for any value of u. For this reason we cannot con- 
sider the self-covariance here. In Part c of this Section we consider a possible 
way of removing the divergence at 0, and consequently, then, we are able 
to calculate the self-covariance for non-vanishing values of u. 

c. Refinement oj the Theory oj Pari b. Here we propose a simple modifi- 
cation of the model of Part b, removing the divergence of the integral of 
S(oi) at oj = 0. The modification considered here, although it is one of several 
possibilities any one of which is sufficient for removing the divergence 
(See Section 6.), is perhaps the only one that is sufficiently simple to treat 
in a memorandum of this scope. The results of this section are thus intended 
to be only provisional and suggestive. 

Let us reconsider the statistics of the function h(r) giving the separation 
between the surfaces near a point r on the mid-plane. The distribution of 
h's considered in the last section is open to several criticisms: (1) it possesses 
no characteristic length parallel to the mid-plane; and (2) the self-covariance 
fh(r)lii(r') (e) does not exist for any value of r — r'. 

To correct partially for these difficulties we replace Eq. (5.22) by 

^"^Mrfwy (5 - 25) 

w 



where ( is a new characteristic length. The self-covariance h\(r)lh(r') 
based upon (5.25) now exists for all values of r — r' except 0. Thus we still 
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have the objection that the variance//? is infinite; however, this will 
cause us no trouble. 
With Eq. (5.25) instead of (5.22) the spectral of density Ri takes the form 

I k k dk 

„ i + l 2 & ' co 2 + D*k* 

= ( x /4*s")-f3a &Arl/c»-Q(y), 



S(fi») = x D/Ts"-0a 1 R i A 2 - f - 



Q(y) = 



y [ y - - log y 



y = fw/D. 



(5.26) 



1 + f 

In obtaining the above equation we have made the usual assumption that 
the frequencies of interest are larger than coo = 4ir 2 D/L 2 , and have replaced 
the original sum by an integral. The function Q(y) has the following proper- 
ties : 



Q(y) ~ — - y log y for y <£ 1 



(5.27) 



Q(y) ~ 1 for y » 1 

Hence for w « D/t 2 , S(u) <* log w, the integral of which converges as o> — > 0; 
whereas, for w » D/F, 5(w) differs negligibly from that given by the unre- 
fined theory (Eq. (5.24)). 

The self-covariance Ri(l) Ri(l + u) now exists for all non- vanishing u and 
is given by 

»oo pi —Duk^kdk 

= ( x /4«") -/?a 2 #4, • e Dul(2 [- Ei( - Du/t 2 )} J 



(5.28) 



where 



-Ei(-x) = J e v dv/v, 

en — log yx for x<Kl, 



~ — for x » 1, 

a; 



7 = 0.5772. 



Thus for u « P /D, Ri(t)Ri(l + u) a - log (7ZVO and for « » ^/Z), 
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Thus we have illustrated how one modification of the model has removed 
the divergence at o> = 0. 

It appears from the treatment here and in part b that roughness and 
diffusion in two dimensions are essential (at least in a linear treatment) 
features in obtaining S(io)a 1/oj. In the case of a non-linear coupling (to be 
considered in a later paper) a "self-induced" roughness effect may occur 
without introducing roughness ab initio as an intrinsic feature of the model. 

6. Summary 

(a) If the resistance deviation R x (t) is related to the concentration devia- 
tion Ci(r, /) of a diffusing medium (particles or heat) by the linear functional 



RiO) = / f(r)cdr, t) dr, (6.1] 



where r is a vector and dr a volume element in a v-dimensional space of 
volume A, , then the spectral density 5"(o>) of Ri(t) is 

v 2 X AJ) k 2 \f k \ 2 
ir s" or + D L ft* 

where D is the diffusion constant, s" is defined by Eq. (3.14), x is the Boltz- 
mann constant, a is the frequency (in radians per sec), k is the wave number 
vector in ^-dimensional ft-spacc limited to a discrete lattice of points (denned 
by Eq. (3.2)) over which the summation is taken, and/* is the ftth Fourier 
component of /(r) (Eq. (3.3)). 

(a) If the important terms in (6.2) vary slowly between lattice points in 
/z-space (true if o > w given by Eq. (3.18)), then (6.1) can be replaced by 
the integral 

«, , r +i ,-i xD f \f(k) fVdk ,, ,. 

where the integration extends over the entire /e-space and where f(k) is 
given by (Eq. 3.20) 

f {k) = rTxJ f(r)e- ihr dr. (6.4) 

(b) Let w' be a frequency in the middle of a wide range. Suppose | f(k) | 2 
averaged over the total solid angle in v-dimensional ft-space is proportional 
to k ", where n is an integer, in a wide range of k with k = -y/u'/D in its 
middle. It follows then that S(u) « D~"~"' 2 w - 1+n+, ' /2 a s long as -1 < 2w 
+ v + 1 < 3. As a consequence, we see that with n an integer (as is true for 
the simple cases considered in Section 5) v must be 2 — the only even di- 
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mensionality— in order that -S(w) be inversely proportional to o> in agree- 
ment with experiment. In this case the only allowed value of n is — 1. 

(c) From (b) we have the interesting result that 5(co) is independent of D 
when it is inversely proportional to co. This means that very slowly diffusing 
substances can contribute as much to contact noise as rapidly diffusing 
substances. This result can be derived on quite dimensional grounds and is 
not dependent upon the special assumptions underlying our treatment. 

(d) A system comprising a high resistance layer modulated by the three- 
dimensional diffusion of particles or heat gives S(o)) « or m . See Case a.(t) 
in Section 5. 

(e) In a system composed of a localized contact disturbed by a diffusing 
surface layer (See Case a.(ii), Section 5), the self-covariance Ri(t)Ri(t + u) 
is inversely proportional to A + Du where A may be considered the contact 
area. We have S(a>) « — log a + const, for to « D/A and S(o)) cc co -2 
for to » Z)/A. 

(f) In a system involving the contact between relatively large areas of 
rough surfaces covered with diffusing surface layers (Cases b. and c, Section 
5), we have been successful in obtaining S(a) « or 1 , and also in obtaining 
a reasonable dependence upon the average resistance. 

APPENDIX I 

Spectral Density and the Self-Covariance 

Here we consider in detail the spectral density, the self-covariance, and 
the relation between these two quantities, first for the case of a single random 
variable. The treatment is subsequently extended to the case of a set of 
random variables which necessitates the consideration of the spectral density 
matrix and the covariance matrix. 

Let y(/)be a real random variable whose time average vanishes, y(t) = 0. 
Now the m.s. value of y can be defined 

y\t) = L : .m - [ y\t, t) dt (1-1) 

T-»0O T J— oo 

where y(t, t) = y(t) in the interval — - < t < - and vanishes outside this 
interval. Evidently y((, t) can be expressed by the Fourier integral 

y(l, t) = [ s(«, T)e ial du (1-2) 

•'-00 

where 

1 r +<T/2) 

*(«, r) - — y(t)e- ial dt. 

L-K J-( r /2) 



LINEAR THEORY OF FLUCTUATIONS 139 

By Parseval's theorem we obtain 

-+CO -+0O 

/ y 2 (t, t) dl = 2v J | y(u, t) | 2 da, 



which, when combined with (1-1), gives finally the desired result (using 
the fact that | y(u, t) | 2 is an even function of co) 

Til) = I F(«) da, (1-3) 

Jo 

where 

F(«) = 4tt Lim - | y(co, t) f (1-4) 

T-»0O T 

is the spectral density. 
By a procedure not very different from the preceding, one can show that 

y(l)y(t -\- u) = / F(co) cos cou do>, (1-5) 

Y(oi) = — / y(t)y(l + u) cos uu du. (1-6) 

■k J 



The quantity (y(l)y(l + u) is called the self-covariance. 

Now let us suppose that we have a set of random variables y,(/) which are 
in general complex and whose time averages vanish. We are then led to 
consider, instead of (1-3), relations of the form 

y.(/)y*(/) = f YM da (1-7) 

Jo 

where now 

Yij(co) = 2ir Lim - [y,-(«, r)y*(w, r) + y,(— w, r)y*( — w, r)] (1-8) 

T-«oo T 

in which 

y.(«, t) = — / yi(l)e tat dt. 

ZlT J-(t/2) 



Instead of self -co variances like y(l)y(l + «) we have to consider a covari- 

ance matrix of the form y ,■(/))',■ + *<)■ Since we shall not have occasion in 
this paper to consider the relation between the spectral density matrix and 
the covariance matrix we will not consider the derivation of the analogue of 
Eq.(I-5). 
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APPENDIX II 
Thermodynamic Theory of Fluctuations 



The value of the quantity Ci(r, t)c t {r\ t) or (<*(*)<£»(')) is determined from 
equilibrium considerations. Before going into the above continuum problem 
let us first consider the problem for the case of a system described by a finite 
set of variables. More specifically let us suppose that the state of the system 
subject to certain restraints (i.e. fixed total mass and energy) is described 
by the set of variables %i , ■ ■ • , x„ . Let the equilibrium state be given by the 
values xi , • • • , x„ , and let 

Xi = *S + «. . (II-l) 

If the system is constrained to constant average energy E, the entropy of the 
non-equilibrium state S = S° + AS will be less than S°, the entropy of the 
equilibrium state, by an amount 

AS = -iZ&f <**«», (n-2) 



where 



W&9/V-2 T° \dXidXj). 



Obviously, AS must be the negative of a positive definite quadratic form, 
otherwise the equilibrium state would not be a state of maximum entropy. 
The probability distribution 26 for the a's is given by 

P(a 1} -,«„) = Ne As,x (II-3) 

where N is a normalization factor and x is the Boltzmann constant. Averag- 
ing the products on atj we find that 

/_, Saaiak = X^.t- (II-4) 

i 

Multiplying (II-4) by the arbitrary set ji and summing over i we get 

zlyiSija,ak = xik. (H-5) 

The generalization to a system described by a continuous set of variables 
is not difficult on the basis of (II-5). Now suppose that, in a v-dimensional 
space Ay , we have a system whose state at time / is defined by the continuous 
set of values of the variable c(r,l) = c + ci(r, I) ; we have 

AS = -I [ s"c\{r,t)dr (II-6) 

» H. B. G. Casimir, Rev. Mod. Phys. 17, Nos. 1 and 3, 343-4 (1945). 
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where 



,dc- /c=3 






when s and e are the entropy and energy, respectively, per unit volume (of 
j'-dimensional space). In calculating (II-6) it was assumed that 

f ci(r, /) dr = 0, 

J A, 

expressing the fact that the system is closed. In order to put (II-6) into a 
form strictly analogous to (II-2) we write it 

AS = -i f f s"5(r - r') Cl (r, l) Cl (r', l) dr dr'. (II-7) 

2 J a, Ja, 
We see that the equation analogous to Eq. (II-5) must be 

( [ y{r')s"8(r - r"Mr", t) Cl (r, t) dr dr' = X y(r) (II-8) 

Ja, Ja, 

where y(r) is an arbitrary function. Integrating (II-8) with respect to r' f 
and using the fact that the delta function is denned by 

Jy(r')8(r' -r) dr' = y(r) 

we readily arrive at the result 



ci(r, t) Cl (r', t) = i 8(r - r'). (II-9) 

Using the Fourier space-expansions of C\ and 5(r) 

ci(r, t) = E' c k (l)e ikr , 
k 

Ay k 

in the region A v = L X X ■ ■ ■ XL, with &, = 2inii/Li , we can write (II-9) 
over into the equivalent expression: 



c k {t)c k {l) = -^-rJkk', (1M0) 



where 

fl if k = k', 
8kk' = \ 

otherwise. 



